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Abstract. Spectral properties of a coupled N x N potential model obtained with the 
help of a single non-conservative supersymmetric (SUSY) transformation starting from 
a system of N radial Schrodinger equations with the zero potential and finite threshold 
differences between the channels are studied. The structure of the system of polynomial 
equations which determine the zeros of the Jost-matrix determinant is analyzed. In 
^ particular, we show that the Jost-matrix determinant has N2'^~^ zeros which may all 

■ correspond to virtual states. The number of bound states satisfies < nt, < N. The 

, maximal number of resonances is — {N — 1)2^~^. A perturbation technique for a small 

coupling approximation is developed. A detailed study of the inverse spectral problem is 
given for the 2x2 case. 

(N 
O 

a^ 

O ■ PACS numbers: 03.65.Nk, 24.10.Eq 

^ ' 1. Introduction 



Almost all low-energy collisions of microparticles with an internal structure (i.e., atom- 
atom, nucleus-nucleus etc) include inelastic processes such as excitations of internal degrees 
of freedom of colliding particles or processes with rearrangements of their constituent parts. 
These processes can be described by a matrix (more precisely multichannel) Schrodinger 
equation with a local matrix potential [1, 2]. One may be interested in both direct and 
inverse scattering problems for this equation. The method of SUSY transformations is 
known as a powerful tool for solving both types of problems for a single- channel Schrodinger 
equation [3]. Nowadays, the first attempt to generalize the method for a coupled- channel 
Schrodinger equation with different thresholds is given in [4, 5]. This attempt is based on 
a non-conservative SUSY transformation (contrary to [6, 7]), i.e. a SUSY transformation 
that does not preserve a boundary behavior of solutions. The main advantage of such 
transformations is a possibility to obtain multichannel potentials with a non-trivial coupling 
starting from the zero potential. 
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The present work is aimed at the investigation of spectral properties of these SUSY 
potentials. Our approach is based on an analysis of the Jost matrix. In the non-relativistic 
scattering theory the Jost matrix plays a fundamental role similar to the scattering matrix. 
The zeros of the Jost-matrix determinant define positions of the bound/ virtual states and 
resonances [1,2]. Therefore, studying the zeros of the Jost-matrix determinant allows one to 
analyze the spectrum of the model. A closed analytical expression of the Jost matrix, as well 
as potential, resulting from a non-conservative SUSY transformation of the zero potential 
is obtained in [4]. The analysis of spectral properties for such potentials was not presented 
up to now despite the fact that the Jost matrix is well known [8]. This may be explained by 
the fact that the spectrum of the potential after a non-conservative SUSY transformation 
changes essentially and to find these changes one has to find all the zeros of the Jost-matrix 
determinant. More precisely, no one spectral point of the initial Hamiltonian belongs to the 
spectrum of the transformed Hamiltonian. As a result, a supersymmetry algebra, which 
is always present in the case of conservative SUSY transformations, cannot actually be 
constructed here and the word 'SUSY transformation' is only a formal heritage from the 
previous conservative case [6, 7]. 

The principal point of this paper is to show that the qualitative behavior of the 
spectrum of (non-conservative) SUSY partners of the vanishing multichannel potential 
with threshold differences may be studied for an arbitrary number of channels, A^. We 
think this is a very strong result, since even for the case N = 2 the full analysis of the 
spectrum is a very complicated problem [8, 10, 11]. The main reason for this is an extremely 
rapid growth of the order of an algebraic equation defining the spectrum with the growth 
of the number of channels. 

The paper is organized as follows. We start with preliminaries, where we give basic 
definitions and equations. Section 3 is devoted to the analysis of the number of bound 
states resulting from a non- conservative SUSY transformation of the zero potential as a 
function of the parameters defining the transformation. This analysis is based on the study 
of the properties of the eigenvalues of the Jost matrix. Following similar lines we analyze 
the possible number of virtual states in section 4. Once the bound and virtual states are 
analyzed we can formulate conditions under which resonances may appear; this is made 
in section 5. The behavior of the Jost-matrix determinant zeros is studied in section 6 in 
the approximation of a weak coupling between channels. In section 7 we deal with the 
particular two-channel case. In this case we express parameters of the potential in terms 
of zeros of the Jost-matrix determinant, i.e. solve an inverse spectral problem. The main 
results are summarized in the conclusion. 

2. Preliminaries 

Let us first summarize the notations used below for coupled-channel scattering theory 
[1, 2, 9]. We consider a system of coupled radial Schrodinger equations for the s- waves 
that in reduced units reads 

H'ijj(k,r) ^ K^'ijj{k,r), r e (0, oo) (1) 
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with 

H = -l^^^V(r). (2) 

where r is the radial coordinate, V{r) is a,n N x N real symmetric matrix, 1 is the unit 
matrix, and ip may be either a matrix- valued or a vector- valued solution. By k we denote 
a point in the space C^, k — {ki, . . . , k^}, /cj e C. A diagonal matrix with non- vanishing 
entries ki is written as K — diag(A;) = diag(/ci, . . . , fcjv)- The complex wave numbers k^ are 
related to the center-of-mass energy E and the channel thresholds Ai, . . . , A^v, which are 
supposed to be different from each other, 7^ Aj, by 

k]^E-Aj, Ai = 0. (3) 

We assume here that Ai = and the different channels have equal reduced masses, a case 
to which the general situation can always be formally reduced [2] . 

Let us recall basic definitions from SUSY quantum mechanics [3, 4, 5, 6, 7]. It is 
known that the solutions of the initial Schrodinger equation (1) may be mapped into the 
solutions of the transformed equation with help of the differential-matrix operator 



ip{k, r) — Lip{k, r) 



4,(k,r). (4) 



The transformed Schrodinger equation has form (1) with a new potential 

V(r) = V(r) - 2U'(r) . (5) 

Matrix U is called superpotential 

U{r) = rj'{r)r]~^{r) , (6) 

and expressed in terms of a matrix solution rj of the initial Schrodinger equation 

Hr){r) = -lCr]{r) , (7) 

where K, — diag(«;) = diag(«;i, . . . , ktv) is a diagonal matrix called the factorization 
wave number, which corresponds to an energy £ lying below all thresholds, called the 
factorization energy. The entries of /C, thus, satisfy £ = —nf + A^; by convention, we 
choose them positive: > 0. Solution rj is called the factorization solution. 
In the case of the zero potential V" = 0, 77 contains only exponentials 

r}{r) = cosh(/Cr) + /C"^ sinh(/Cr)C/o . (8) 

The symmetric matrix Uq is the superpotential at r = 0, which can be chosen arbitrary. It 
is convenient to introduce special notations aj for the diagonal and for the off-diagonal 
entries of f/o. 

The Jost matrix of a (non-conservative) SUSY partner of the A^-channel zero potential 
reads [4] 

F{k) ^ {}C-iK)-\Uo-iK) , (9) 
which is also the Jost matrix obtained in [8]. 
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The necessary and sufficient condition on the parameters (factorization energy S and 
superpotential at the origin Uq) to get a potential without singularity at finite distances is 
obtained in [10, 11]. This condition is the positive definiteness of matrix /C + Uq: 

/C + t/o > , (10) 

which puts some upper hmit on the factorization energy S at fixed Uq. 

Zeros of the Jost-matrix determinant define positions of the bound/virtual states and 
the resonances. Thus, to find these positions we have to solve the following equation 

detF{k)^0, (11) 

taking into account the threshold conditions (3). According to (9), the roots of equation 
(11) are defined by the roots of 

detS(A;) = 0, Kj-ikjy^O, j^l,...,N, (12) 

where 

B{k) = Uo- iK . (13) 

In what follows we concentrate on the analysis of the zeros of det-B only keeping in mind 
that some of them may be cancelled in detF if kj = —inj. Our starting point is thus a 
system of algebraic equations (12) and (3) which reads, with certain coefficients a], 

N N N 

{-ifkik2 ...kN + Y^ a^^_i J]^ fc; + . . . + ^ a{kj + Oq = , (14) 

j=i i=i,ii^j j=i 
k] - A;^ + A,- = . (15) 

First we show that system (14), (15) can be reduced to an algebraic equation of the 
N2^~'^ degree with respect to one momentum, say /ci, only. Indeed, any momentum enters 
equation (14) only linearly. Therefore it can be rewritten in the form 

kNPi{ki, kjsr-i) = Qi{ki, kN-i) , (16) 

where Pi{ki, . . . , /cat-i) and Qi{ki, . . . , /cat-i) are polynomials of the first degree in each of 
the variables ki, . . . , fcjv-i- It is important to note that given all momenta ki, . . . , k^-i this 
equation defines /cjv in a unique way if Pi does not vanish. On the other hand we can square 
the left- and right-hand sides of (16) thus obtaining an equation where kN enters only in 
the second degree and polynomials and Ql are polynomials of the second degree with 
respect to their variables. But in the equation thus obtained using threshold condition (15) 
we can replace all second powers of the variables kj, j = 2, . . . , N hj kf — Aj, which makes 
disappear both variable kj^ and the second power of kj, j = 2, N —1 from the resulting 
equation and raises the power of ki till 2N. We thus sec that after these manipulations 
variable k^^i enters in the resulting equation only in the first degree and the equation can 
be rewritten in form (16) 

kN-iP2{ki, • • • , ^iv-2) = Q2iki, kN-2) , (17) 

where P2{ki, . . . , kN-2) and Q2{ki, . . . , ki^_2) are polynomials of the first degree in each of 
the variables k2, ■ ■ ■ , kN-2- From (17), given ki, . . . , kN-2, not a zero of Pn-2, we obtain 
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/cAr_i in a unique way. We note that the system (17), (16) and (15) where from (15) the 
last equation klf — kf + Aj^ = should be excluded, is equivalent to the original system 
(14), (15). 

It is clear that we can repeat the above process N — 3 times more to get an equation 

k2PN-i{ki) = QN-i{ki) (18) 

and finally 

PN{ki) = (19) 

with Pjv of order N2^~'^. Note, that the subscript in and Qk indicates nothing but the 
step in this procedure. It is evident that any ki which (together with k^, ■ ■ ■ , /cat) solves 
the system (14), (15) is a root of (19). The converse is also true. Indeed, given a root ki 
of (19), but not a root of Pn-i, we find from (18) a unique /c2. Once we know ki and k-z 
wc find A;3 from equation previous to (18) and so on till kj^ which is found from (16). It 
is also clear that in this way we can get N2'^~^ number of sets ki, . . . ,kN (some of them 
may coincide) each of which solves the system (14), (15) so that the same number N2'^^'^ 
is the number of possible solutions of this system and the system (19), (18), . . . , (16) is 
equivalent to the initial system (14), (15). 



3. Number of bound states 



In the following, except for sections 6 and 7, wc will consider all quantities as functions 
of the momentum ki. Other momenta arc expressed in terms of ki from the threshold 
conditions (3). Since in this section we are interested in the number of bound states we 
will consider only the negative energy semi-axis E e (— oo, 0). It happens to be useful to 
change variables kj in favor of kj as kj — ikj and rewrite the threshold conditions (15) 
accordingly 

k, = ^Jkf + AJ, (20) 

where we have chosen only the positive value of the square root since in this section 
we analyze only the point spectrum of H, which restricts all momenta kj to be purely 
imaginary with a positive imaginary part so that kj = \kj\. 

Prom (9) it is clear that all the zeros of detF are at the same time the zeros of the 
determinant of matrix B (13) and vice versa. This follows from (12) and the positive 
definiteness of matrix /C — iK in the momenta region we consider so that neither of the 
roots of det B solves the equation det(/C — iK) — 0. 

Since det B — Y[j=i "^j where Xj are the eigenvalues of B, 

B{h)xj{h)^Xj{h)xj{h), (21) 

the equation deti?(A;i) = is equivalent to Xj{ki) = 0, j = 1,...,A^. Matrix B is 
symmetric with real entries in the momenta region we consider, B — Uq + K — B^, 
which implies the reality of both Xj{ki) and Xj{ki). Here we introduced a diagonal matrix 
K — \K\ — diag(^i, . . . , ^;v)- 
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Another property of Aj(A;i) important for the analysis is their monotony as functions 
of ki that we prove below. 

For a fixed K let us consider a deviation of \j(ki) for a small increment of argument 
^1, i.e. Xj{ki + Ski) — Xj{ki) + SXj{ki) assuming SK — diag(5A;i, . . . , Sk^) real, positive 
definite (since 5kj > 0, Vj) and infinitesimal. Prom (21) one gets 

B{ki + Ski)xj{ki + Ski) = Xjih + 6ki)xj{ki + 6ki) . (22) 

Here according to (13) B{ki + 6ki) = Uq + K + 5K and the increment of B{ki) is just 
5B — SK which plays the role of a small perturbation of B(ki). Therefore we may 
calculate the shifting of the eigenvalues produced by such a perturbation using a (Rayleigh- 
Schrodinger) perturbation theory. Thus, for a non-degenerate eigenvalue Xj the first order 
correction reads 

SXj = {xj\SB\xj) > (23) 

where the inequality follows from the positive definiteness of SB — SK, which in turn 
implies monotony of the eigenvalues as functions of the momenta ^i. For a degenerate 
eigenvalue corrections are obtained by diagonalizing the same perturbation operator SB 
restricted to a linear span of unperturbed eigenvectors corresponding to a given eigenvalue, 
which still leads to positive corrections because of positive definiteness of SB. 

From here it follows that any eigenvalue Xj{ki) may vanish i.e. change its sign, only 
once. Moreover, Xj kj > as ki — > oo. Hence, the number of negative eigenvalues of B 
at ki — 0, i.e. at the energy of the lowest threshold, is just the number of bound states. 
Thus, to count the number of the bound states, Ub, one has to consider the eigenvalues 
Xj(ki), j — 1, . . . ,N ol matrix B{ki) at = 0, 

S(0) = C/o - i diag ii^) = C/o + dmg{^) (24) 

so that 

n,^\{N-K), A = f]A,, A.-^- (25) 

To clarify this formula we notice that in the absence of bound states all Aj = 1, A = 
so that rih = 0. Every bound state is responsible for the change of the sign of only one 
eigenvalue from positive to negative thus raising —A by 2 units, i.e. —A — > —A + 2 with 
nb+ 1. This justifies the factor 1/2 in (25). 
Summarizing, we see that the number of bound states is bounded hy < Ub < N. 
Figure 1 shows the eigenvalues of matrix B as functions of ki for the case A^ = 3. Two 
eigenvalues cross the axis which corresponds to the case of rib = 2. The last comment in 
this section is devoted to equation (10). Now it can be seen that the factorization energy 
should be chosen lower than the ground-state energy for the transformed potential, S < Eg, 
if any. 

4. Number of virtual states 

According to the definition of a virtual state [1, 2] in this section we will need to consider 
the channel wave numbers kj lying both in the positive and the negative imaginary semi- 
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Figure 1. Typical behavior of S-matrix eigenvalues, = 3. The case of two 
bound states with energies Ei = —51.8611 and E2 = —8.8852 is presented. The 
black squares show positions of these bound states. The corresponding parameters 
are oi = -3, 02 = -8, 03 = -1, P12 = 0.5, /3i3 = 0.4, P23 = 1, A2 = 15, A3 = 25. 



axes of the corresponding momenta complex planes and consider the full imaginary axis 
for ki, i.e. ki G (—00, 00). The other momenta, k2, ■ ■ ■ , kj^, belong to either the positive 
or to negative parts of their imaginary axes in agreement with the threshold conditions 

h = ^j\/kl + ^j, ^j = ±, J = 2,...,N. (26) 

Since in (26) all combinations of signs are now possible it is convenient to introduce 
special notations for these combinations. Denote a = (+, ±, . . . , ±) a string of signs 
with cTj being its j-th entry, which corresponds to the sign in (26) for the j-th momentum 
for J > 1. The first symbol "+" in a indicates that all momenta kj are expressed in 
terms of ki. Let n+((T) + 1 and n-{a) be the numbers of "+" and "— " signs in this 
string. We notice the following evident combinatoric properties of n_(cr) and n_^{a). First, 
n_|_((T) + n_(o") + 1 = which implies 

J2 IM^) + n-i^) + 1] = A^2^'' • (27) 

Here and in what follows the summation over a includes all 2^^^ possible sign combinations. 
Next, a symmetry between "+" and "— " leads to the following relation 

5^n_(a) = $^n+(a) = (iV - 1)2^"^ (28) 

According to (13) every sign combination leads to its own B matrix defined by the 
corresponding K matrix so that both K and B should carry an additional information 
about this combination. Therefore 

B'' = Uo + K% i?" = diag(fci, cTsfca, • • • , ct^v^tv) (29) 

and we denote Aj(fci), j = 1, . . . , the eigenvalues of B'' . 

In order to find the zeros of the Jost-matrix determinant corresponding to the virtual 
states we should find the purely real solutions of the equations A°"(A;i) = 0, j = 1, . . . , 
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ki ki 

Figure 2. Typical behavior of the eigenvalues Xj{ki), = 3, is shown. 
Each plane corresponds to a particular choice of string a: (a) a = {+ + +), 

(b) (T = (+ H — ), (c) a = (H ), (d) o" = (H Stars, squares 

and circles correspond to the virtual states. Virtual states are denoted by 
the identical symbol if they belong to the same eigenvalue Xj{ki). The 
corresponding parameters are ai = 3, a2 = 5, as = 9, (3i2 = 0.5, /?i3 = 0.4, 
/323 = 0.2, A2 = 15, As = 35. 

for all 2^^^ matrices B"^. Although the k/s are real, but bearing in mind our replacement 
kj = ikj, throughout the text we call these zeros purely imaginary. Finally we note that 
since matrix JC — iK in (9) is not positive definite for an arbitrary a anymore, in some 
particular cases some of the zeros of B may be cancelled by the zeros of det(/C — iK) and 
will not correspond to virtual states. Nevertheless, omitting these particular cases, we will 
concentrate on an analysis of the zeros of det-B only. 

Eigenvalues Xj{ki) are monotonous functions of ki in two cases only: (i) a = 
(+, +,...,+) and ki > 0; (ii) a = (+, —,...,—) and ki < 0. In general, an eigenvalue 
Xj{ki) may have minima/maxima for ki ^ which may lead to two or even more roots in 
equation A°'(fci) = 0. We illustrate this behavior for = 3 in figure 2. The monotonous 
lines in the right/left part of figure 2(a)/(c). correspond to case (i)/(ii). The position of 
the zeros of the eigenvalues is shown by stars, squares and circles. It is clearly seen that 
the total number of the roots of all equations Xj{ki) = is (A^2^~^)|jv=s = 12 which all 
correspond to virtual states. 

A change of parameters may result in shifting the position of the virtual states only 
without changing the number of zeros (i.e. virtual states). For instance, in the simplest 
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case we may shift all diagonal entries of Uq by a number Aq, Uq ^ Uq + XqI, thus shifting 
all eigenvalues of B by the same number, Xj{ki) — > Aj(fci) + Aq. 

Let us consider a specific eigenvalue defined by a string ctq, with a local maximum at 
ki — ki^max-i ^'j" (ki^rnax) = '^j,max- One Can always shift all the eigenvalues by the value 
^j,max such that the curve AJ°(^i) touches the ki axis at the point = k\,max nieaning 
that ki^max not only becomes a root of the equation AJ°(^i) = but this root is multiple 
(of multiplicity 2) and by a small additional change of other parameters it can be split into 
two simple but complex roots. This is just in this way two virtual states collapse producing 
a resonance; a subject which deserves a special discussion (see the next section). Pairs of 
virtual states which may collapse are shown in figure 2 by squares and circles. 

It is not difficult to convince oneself that for any given f3ji the situation when all the 
zeros of the Jost-matrix determinant are purely imaginary may be realized by a proper 
choice of aj. To see that let us consider the asymptotic behavior of AJ for — > oo, when 
all off-diagonal entries of B become neghgibly small, 

A^ ~ ^1 + ai , (30) 
AJ ~ aj^kl + Aj + aj = aj (^\h\ + ^ + . . + . (31) 

l^il^oo. (32) 

Numbers n+{a) and n-{a) determine the corresponding numbers of increasing and 
decreasing eigenvalues at positive infinity. The eigenvalue AJ increases both at negative 
and positive infinity. Now if we choose all aj sufficiently large in absolute values and 
negative we can always guaranty the location of a root of the equation A^(A;i) = near the 
point ^1 = «! and at the same time the location of two roots of the equation X^{ki) — 
with corresponding aj = + near the points ki = ia^. Thus, for each cr we can obtain 
2n+((T) + 1 zeros. The total number Uy of these zeros may be calculated by formulas (27) 
and (28) 

Uy = J][2n+(c7) + 1] = 7V2^-i , (33) 

which coincides with the total number of all possible roots of the system (14), (15) and is 
just the maximal possible number of virtual states. Hence, in this case all the roots are 
purely imaginary. In the next section we consider the case when some of the zeros may 
merge, become complex and produce resonances. 

5. Number of resonances 

For simplicity independently on whether or not it can be seen in a scattering we call any 
pair of complex zeros k — ±kr + iki oi the Jost-matrix determinant a resonance keeping 
in mind that to be really visible in a scattering a resonance behavior of the corresponding 
cross-section should be narrow and sharp enough. 

Conservation of the number of zeros of an n-th order algebraic equation under a 
variation of parameters included in its coefficients, which keeps unchanged its order (in our 
case this is equation (19) obtained from the system (14), (15)) applied to our case leads to 



Non-conservative multichannel SUSY partners of the zero potential 



10 



the following relation rih + -\- 2nr = N2 , where Ub, Uy and Ur are number of bound 
states, virtual states and resonances respectively. The aim of this section is to establish 
the maximal number of possible resonances accepted by a non-conservative SUSY-partner 
of the vanishing potential. 

Evidently, the maximal number of resonances corresponds to the minimal number 
of bound Ub and virtual states. These numbers would both become zero if no one of 
the B matrix eigenvalues intersected the ki axis. But as it was noticed in the previous 
section there always exists an eigenvalue with the asymptotic behavior given in (30), 
i.e. ranging from — oo to +00 and, hence, it intersects ki axis always and for all possible 
values of a. We thus see that the minimal number of real zeros that all eigenvalues may 
take is achieved if all eigenvalues Aj(A;i), j > I are nodeless and curves Xi{ki) intersect 
^1 axis only once for every given sign combination a. To realize this case, we should 
choose parameters included in Uq in a such way that the global minimum Xj^rnin every 
eigenvalue AJ(^i) with Uj — -\- (they tend to -|-oo when |^i| — > 00) be positive AJ_^j„ > 
and, respectively, global maximum Xj^^ax every eigenvalue AJ(^i) with aj — — (they tend 
to —00 when |^i| — > 00) be negative X'j^^nax < 0- Under these conditions only eigenvalues 
A^(A;i) have zeros. The possibility that these eigenvalues have only one zero can always 
be realized. This can be demonstrated for small enough values of jSij (so called weak 
coupling approximation, see the next section) which in the limit Pij = for all i,j gives 
a very simple behavior of the eigenvalues. For instance, for A^+i — Aj large enough and 

min(A/AT + a,) > max(— a/AT + a,) the straight line A'i'(A;i) never intersects with the 

j _ j 

hyperbolas Aj(A;i) so that small perturbations coming from small non-zero /3j7-values (in 
a physical terminology these perturbations shift the zero width resonances from the real 
energy axis to the complex plane) do not change the monotonous behavior of Xi{ki) and, 
hence, do not bring additional roots to the equation A^ (^1) = 0. 

Thus, we see that the minimal value of virtual states with the absence of bound states 
is equal to all possible sign combinations of a which is nymin — — 2^~^. Hence, the 

maximal possible number of resonances is obtained by subtracting this number from the 
total number of solutions, i.e. 

2nr,max = N2^-' - 2^'^ = (TV - 1)2^-1 . (34) 
6. Weak coupling approximation 

For the number of channels N > 2 there is no way to get analytical solutions of system 
(14), (15) but if the coupling parameters Pij are small enough assuming analyticity of 
the roots of the Jost-matrix determinant as functions of Pij a perturbation technique may 
be developed. In this section we demonstrate this possibility by obtaining first order 
corrections to unperturbed values of the roots of the Jost-matrix determinant corresponding 
to Pij = 0. 

For the zero coupling, matrix Uq becomes diagonal Uq — diag(Q;i, a2, ■ ■ ■ , ol^) and the 
system (11), (3) reduces to 

(«! - iko^i){a2 - iko,2) ■ ■ ■ (oiiv - iko,N) = , (35) 
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/c2,-/co'i + A, =0, j = 2,...,N, (36) 

where the additional subscript corresponds to the uncoupled case. Its solutions have the 
form 

^^I/"^^ = - , ko^m^ = (Jrn\J -a% + A^v - A^ , TTl^ N , (37) 

where m = 1, . . . ,N. Let us explicitly indicate the meaning of sub- and superscripts in 
(37): the second subscript m in k^^'^'' corresponds to the channel, the first superscript 
j indicates a row number in (37) and a indicates one of all 2^~^ combinations of signs. 
Thus, we see once again that the total number of solutions of the system is N2^~^ and it 
does not depend on whether or not the coupling is absent. Note that every energy level 
Ej — —a'^ + Aj corresponding to a row in (37) is 2^~^ fold degenerate. Below we show that 
under a small coupling every degenerate level Ej splits by 2^""^ sub-levels and we will find 
approximate values of the splitting. But the unperturbed j'-th momentum corresponding 
to this level simply equals kQ'^^ = —iaj. Therefore, instead of our previous convention to 
express all quantities in terms of ki, it is convenient here to express corrections to the j-th 
momentum produced by a perturbation in terms of unperturbed j-th momentum k^^^'f^ ■ 
This is always possible due to the fact that all momenta have equal rights. But now we 
have to change our signs convention introduced in section 4 where the first momentum 
ki entered in the string a always with the positive sign (cri = -|-). Now we have j-th 
momentum kj e (—00,00) and aj = -|- in string a. 

Prom (37) we learn that no coupling implies no finite-width resonances but as we 
discuss below the zeros lying above the first threshold may be associated with zero-width 
resonances which acquire a non-zero width under a small coupling. 

From the first row of (37) we conclude that the corresponding 2^^^ zeros with 
El = —a\ arc always below the first threshold (bound or virtual states). Energy 
En = + A„, n = 2, . . . , N, may be positive with respect to the first threshold and just 
these {N — 1)2-^"^ zeros are associated with the zero-width resonances. According to our 
convention a resonance corresponds to a pair of complex zeros. Here we can easily compute 
the number of the zero-width resonances, Uzwr, which is Uzwr — {N — 1)2^~^ which agrees 
with the maximal number of possible resonances obtained in the previous section. 

Unperturbed B'^ matrix we denote Bq is diagonal 

Bq = diag(Q;i + aiki, a-i + (72^2, . . . , ctiv + o^kj^) (38) 
and its eigenvalues Aq ^ coincide with its diagonal entries 

= % + h = KM = + (Ji^k^j+^i-^j. (39) 
Z = 1,...,7V, l^j. (40) 

For simplicity we assume all coupling parameters /3y proportional to the same small 
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parameter (3 so that the perturbed B" matrix reads 

B'^^B-^+pB, B^Wh^iW, bjj^O, j^l,...,N. (41) 

Now as it was mentioned above assuming analyticity of eigenvalues of this matrix as 
functions of /? we can develop them in a Taylor series with respect to /3, 

- K,j + + + ■■■, (42) 

where the first subscript number is just the power of (3. First we notice that the 
perturbation B has zero diagonal entries which results in A^j = 0. To get the second 
order correction we are using the usual Rayleigh-Schrodinger perturbation approach which 
leads to 

N ^2 

In what follows we also assume that we can neglect the higher-order corrections to the 
eigenvalues. 

Actually, our aim is to find corrections to the unperturbed degenerate j'-th Jost-matrix 
determinant zero given in (37). Assuming a Taylor series expansions for this root over the 
small parameter j3 indicating it now explicitly 

kj = k^'f + pci + p^C2 + ... (44) 

we find coefficients Ci and C2 from the equation 

= KAk) + ^Uh) = Q ■ (45) 

For that we develop X'jikj) in a Taylor series in (3 parameter considering its P dependence 
as given through kj and (44). The term (43) contains the factor /?^, therefore in its 
denominator we simply put ^qJ instead of kj. The Aij-dependence of the term XQj{kj) 
is given by (39) and its /3-dependence is obtained via (44). Thus, the left hand side of 
equation (45) is presented as a series over the powers of (3 where every coefficient should 
vanish. This leads to Ci = and 

C2= E 1 ' (46) 

i=i,ij^j ai + aiJa^j + Ai- Aj 

Finally up to the second order in (3 we obtain the roots of system (37) 



A:?-^ = -^a2 + ^ E fftr^ , k^^-^ ^ aj (k^^-A\ - 



1=1 <^i+<^iy <^N+'-^i~'-^N \ \ / 

Here each row is obtained by applying equations (43), (44), (45) and (46) for j = 1, . . . , A^, 
respectively, and m = 1, . . . , N, m ^ j for each j. The square roots in the last column of 
(47) should be expanded in Taylor series up to (3'^. 
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From here it is easily seen that, when < A^, purely imaginary unperturbed zeros 
km = —io-m move from the axes to the complex plane due to the real part of corrections. 
For instance for /c2, the real part reads ±/3^y'A2 — a^/ {a\ — a\-\- A2). We thus confirmed 
the previous statement that zero width resonances acquire non-zero widths. 



7. Zeros of the Jost-matrix determinant for N — 2 

The particular case of two coupled channels is important both from practical and theoretical 
point of view. Let us recall the following inequalities for the number of the bound/virtual 
states and resonances obtained in sections 3, 4 and 5: < n?, < 2, < < 1, < < 4. 
The same inequalities are obtained ior N — 2 in [10, 11] from another approach. The 
two-channel problem is the only one where one is able to get analytic expressions for the 
Jost-determinant roots, i.e. to solve the direct problem consisting in finding the positions 
of the bound/virtual states and resonances. This possibility is based on the fact that the 
roots of the algebraic equation of fourth, (A^2-^~^)|jv=2 = 4, order may be expressed in 
radicals. Thus we obtain zeros as functions of parameters defining the potential. One may 
be interested in solving the inverse problem: to express parameters of the potential from 
the knowledge about positions of zeros of the Jost-matrix determinant. In principle, one 
may try to inverse radicals, but we propose a more elegant way below. 

To simplify the notations, we choose in this case A2 = A 7^ 0. The potential, which 
is known as the Cox potential [8] , depends on three parameters appearing in matrix 



Uo- ' (48) 




and on the factorization energy S which is upper bounded. The Jost-matrix determinant 
reads 

f{ki, k2) = detF /ci, /ca = — , ^ ■ \ — ■ 

[ki lKi){k2 + IK2) 

The system of equations (14), (15) in this case takes the simplest form 

kl - kl = A, (50) 

{ki + iai){k2 + ia2) + = (51) 

and can be reduced to a fourth order algebraic equation with respect to ki 

kl + iaikf + a2k\ -\- ia^ki + — . (52) 

The coefficients Oj, i = 1, . . . , 4 are given exphcitly in [10], (33a-d). Momentum k2 can be 
found from 

k2{iki — Oil) = Oi2{ki + iai) — , (53) 

which is a direct implication of (51). Equations (52) and (53) are particular case of the 
system (19), (18), . . . , (16) for = 2 in accordance with our general discussion in section 2. 

Let us assume we have found two of the roots of system (50), (51) we denote {k^i \ k^^) 
and {k^\k^^)^ which clearly are functions of parameters ai and 0:2. Their dependence on 
parameters (5 and A is not important for the moment, since both (3 and A assumed to be 
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fixed. Being put back to (51) the equation reduces twice to identity for any values of cti 
and a2-, wliich we write as 

(A;i')+mi)(4')+m2)+/3' = 0, (54) 
{kf^ + ia^){kf^ + ia2) + = . (55) 

The reason why we replaced the identity sign by the equality sign is that these equations 
may be considered as an implicitly written inverted dependence of «i 2 on the set of 
parameters /ci 2 • We may thus fix arbitrary values for and find from (54), (55) 

«! and a2 in terms oik\2 which is a much easier task than finding an explicit dependence 

(1 2) 

of k\'2 

on Q!i and 0:2. For that we have to solve, e.g. for cti, the following second order 
equation 



a\-a,i{k?+kf^)-k'^hfH(5'^ 



0, 



(56) 



with Ri = k\ 
here we find 



(2) 



k^i^ and R2 = kf^ 



k^^ which easily follows from (54) and (55). From 



Q!2 



i{k^^^ + kf^) ± ^f^R[^^4pmjR2 
i{k^'^ + A;f ) T ^-Rl-A(5m2/Ri 



(57) 
(58) 



The upper (resp., lower) sign in (58) corresponds to the upper (resp., lower) sign in (57). 

Q 2) fl 2) 

The values oik\' ' and k^ ' should be chosen so as to warranty the reality of parameters 
a;i,2- 

Once two roots are fixed, (52) reduces to a second-order algebraic equation 
62(^1) = for the two other roots kf^'' and k^^^ thus providing an implicit but rather 
simple mapping between the roots of system (50), (51) and the set of parameters 
(ai,Q;2,/3). Polynomial 62(^1) is the ratio of the polynomial appearing in (52) and 
r^ik,) = kl- ki{kf^ + M'^) + kfhf\ i.e., 

k\ + iaikl + a2kl + ia^^ki + 04 = 7^2(^1) 22(^1) • 

From here we find, with the explicit expression for coefficients Cj, 

Q2(A;i) = (fci + iaif + ki{kf^ + M'^) + 

{2ia, + kf^ + k'tW? + M'^) + «^ - A - k?kf^ 



1,...,4 [10], 



and, hence. 



(3) 



kl 



(4) 



^-Rl-AI3mjR2 + ^i 
Tt\/-Rj-AP^Ri/R2 - Vd'i 



(59) 
(60) 



where Di^ Rl + 4/3^1^ + 4A;f ^/cj^^. The sign before the first square root in (59) and (60) 
should be chosen in accordance with the signs in (57) and (58). 

To find ^2^'^^ we do not need to solve any equation. We simply notice that the equation 
detF{ki, k2) = is invariant under the transformation ki /c2, ai <-> 0:2, A <-> —A. This 
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Table 1. Possible mappings between some experimental data and the Cox potential 
parameters. 



Experimental 


Fixed 


Free 


Restrictions 


data 


parameters 


parameters 






«! , 02 




/3 > y/~krPr 


A ,Eb = ~\l ,Er,Ei 


ai , 0:2 , /9 




Ki > Xb 


A , £"1,2 = -Af^2 






Ki > X2 > Ai 


A,E, = -Xl 


Oi2 




Ki > Xb 



means that being transformed according to these rules equations (59) and (60) give us the 
k2 values: 



1.(3) _ - 

1^0 — "Z" 



1 



1.(4) _ 1 
"2 



4/3^2/ Ri + VD2 



(61) 
(62) 



where L»2 = i^i + 4:p^f- + ^kfh'i\ 



Two initial zeros {kf^ , k^^), (/c}^^ k2^') and threshold difference A are assumed to be 
known from the experiment. For instance, these zeros may correspond to a visible Feshbach 
resonance or two bound states. The possible cases for initial data arc summarized in 
Table 1. The first row of Table 1 corresponds to the case where the position of the 
resonance is known (see section 7.1 below). The second row corresponds to the case where 
the positions of both the resonance and one bound state are known, which allows one to fix 
maximal number of parameters. The third row corresponds to the case where the positions 
of two bound states are known (see section 7.2 below). The last row corresponds to the 
special case when only one zero may be fixed from experimental data. The free parameters 
in Table 1 allow either for isospectral deformations of the potential or for fits of additional 
experimental data as, e.g., scattering lengths (see e.g. [10, 11]). The restriction on the 
factorization energy is deduced from the regularity condition of the potential (10). The 
restriction on the coupling parameter /? is explained below (see (67) in section 7.1). Let us 
now consider examples corresponding to the first three rows of Table 1. 



,(2) l(2) 



7.1. One resonance. 

A resonance corresponds to a pair of complex roots k^^ and k^^^ of the Jost-matrix 
determinant such that ik^^ and ik^^ are mutually complex conjugate. Therefore we 
assume equations (50) and (51) to have two complex roots. Let us define their first-channel 
components as 

fcW =kr + iki , fcf ^ = -kr + iki, fci e M , KeR, K > , (63) 

and write the corresponding energies, (j^^'^^^ , as E^. ± iEi, where we also assume E^ > 
(which means that the upper sign corresponds to k^^ or k^\ depending on the sign of k^). 
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We would like to choose as parameters the threshold difference A, as well as the real and 
imaginary parts of the resonance complex energy, Er,Ei. As exemplified below, these can 
correspond to physical parameters of a visible resonance in some (but not all) cases. In 
terms of these parameters, kr and ki are expressed as 



k-p — . — 



El. 

V2 



E^ + Ef-Er 



-1/2 



E^ + Ef-Er 



1/2 



(64) 



In the second channel the roots 



k 



(1) 



Pr + iPi 



k 



(2) 



-Pr 



can be found from the threshold condition yielding 

1 



Pr 



Pi 



V2 

E^ 
V2 



^Ef + {Er-Ay + Er-A 



nl/2 



n-1/2 



Ef + {Er -Ay + Er-A 



(65) 



(66) 



The upper (resp., lower) sign in (64) corresponds to the upper (resp., lower) sign in (66), 
which means that, for a given zero, the signs of ki and Pi are opposite. Moreover, equations 
(64) and (65) show that, for a given zero, the signs of kr and pr are also opposite. This 
implies that, for the Cox potential, the complex resonance zeros (or scattering- matrix 
poles) are always in opposite quadrants in the complex ki and k2 planes. This has 
important consequences for physical applications: for a resonance to be visible, one of 
the corresponding zero has to lie close to the physical positive-energy region, i.e., close 
to the real positive ki axis and close to the region made of the real positive k2 axis and 
of the positive imaginary k2 interval: [0, iv^]- Consequently, the only possibility for a 
visible resonance with the Cox potential is that of a Feshbach resonance, only visible in 
the channel with lowest threshold, with an energy lying below threshold A. At higher 
resonance energies, the corresponding zero is either close to the fci-plane physical region 
(and far from the /c2-plane one) or close to the A;2-plane physical region (and far from 
the fci-plane one); it cannot be close to both physical regions at the same time, hence it 
cannot have a visible impact on the coupled scattering matrix. Here, we illustrate the 
case of a visible resonance, which is the most interesting from the physical point of view. 
It corresponds to the lower signs in (64) and (66), with a resonance energy E^ such that 
< Er < A, and a resonance width F = 2Ei such that Ei < E^. 

Note, that for non-zero values of the parameters k^ and pj. (which have opposite signs), 
the coupling parameter j3 cannot be infinitesimal: because ai and a2 have to be real, /? is 
restricted to satisfy the inequality 

(3 > V-krPr ■ (67) 

To get a potential with one bound state at energy — A^, we choose the lower signs in 
(57), (58). We then get for k'^\(3) an expression similar to (59), (60), from which the value 
of (3 can be found by solving the bi-squared equation k\ {(3) = i\b- 

Let us now choose explicit parameters. First, we put A = 1. To get a visible resonance, 
we put Ej. = 0.4, Ei = 0.01 (which corresponds to a resonance width F = 0.02), and 




Figure 3. The Cox potential without bound state and with one visible resonance 
of energy = 0.4 and width T = 0.02, for A = 1 and (3 = 0.1 (first row, solid 
lines for Vu and V22 + A, dashed line for V12), with the corresponding partial cross 
section (second row) and phase shifts (third row) for (a) ki = 0.5; (b) ki = 0.7; 
(c) Ki = 1. 



(3 = 0.1. Using (57), (63) and (64), one finds ai = 0.76938 and 02 = —0.766853 (we choose 
the upper signs (57), (58)). The factorization energy, £, is not constrained in this case: 
it just has to be negative. The Cox potential with one resonance and two virtual states 
E^i = -0.560473, E^2 = -0.599544 is shown in the first row of figure 3. 

The diagonal elements of the potentials, Vu and V22 + A, are plotted with solid lines, 
while V12 is plotted with dashed lines. Parameter ki is responsible for the isospectral 
deformation of the potential which results in the behavior of the phase shifts. The second 
row of figure 3 shows the corresponding partial cross sections, where the resonance behavior 
is clearly seen, as well as the evolution of the low-energy cross section, which is related 
to the scattering length. The last row of figure 3 shows the corresponding phase shifts 
for the open channel, where a typical Breit-Wigner behavior (see e.g. Ref. [1]) is seen for 
the resonance, as well as the evolution of the zero-energy phase-shift slope, which is also 
related to the scattering length. 

7.2. Two bound states 

Let us now construct a Cox potential with two bound states, and hence no resonance [10]. 
We choose k^^ = O.li and kf'^ = 1.5i for these bound states and, as in the previous example. 
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Figure 4. The Cox potential (solid lines for Vn and 
V22 + A, dashed line for V12) with two bound states at 
energies Ei = -0.01 and E2 = -2.25, for A = 1, /? = 0.1 
and Ki = 1.51. The left (resp., right) graphic corresponds 
to the upper (resp., lower) signs in (57) and (58). 



we put A = 1 and j3 = 0.1. We thus have ^2 — \/l-01i and ^2 — \/3.25i, which defines 
i?2 in (57), (58). Choosing the upper signs in these equations, we find ai = —0.112649 and 
0^2 = —1.79557, while for the lower signs, we get ai = —1.48735 and 02 = —1.0122. The 
corresponding Cox potentials are shown in figure 4. 

8. Conclusion 

A careful study of spectral properties of non-conservative multichannel SUSY partners of 
the zero potential is given. Our treatment is based on the analysis of the Jost-matrix 
determinant zeros. Generalizing our previous results for the two-channel case [10, 11], we 
have shown that the zeros of the Jost-matrix determinant are the roots of an A^2^~^th- 
order algebraic equation. The number of bound states is restricted by the number of 
channels, < nb < N. The upper bound for the number of resonances is {N — 1)2^~^. 
The generalization is based on the analysis of the behavior of the Jost-matrix eigenvalues. 

In general, an algebraic equation of an order higher than 4 has no solutions in radicals. 
As a consequence, there are no exact analytic solutions of a spectral problem for a non- 
conservatively SUSY transformed Hamiltonian with N > 2. Therefore, the problem of 
finding the approximate solutions appears to be actual. Based on the usual Rayleigh- 
Schrodinger perturbation theory for the eigenvalues of the Jost matrix we develop an 
approximate method for finding the zeros of the Jost-matrix determinant in the case of a 
weak coupling between channels. 

An analytical study of the Jost-determinant zeros is carried out for the two-channel 
case which implies an algebraic equation of the fourth order. A suitable factorization of the 
fourth-order polynomial allows us to develop a procedure which solves the inverse spectral 
problem for this case. The effectiveness of the procedure is illustrated by two examples: a 
potential with one resonance and a potential with two bound states. 
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